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f^) , Abstract: Extended Thermodynamics is the natural framework in which to 

study the physics of fluids, because it leads to symmetric hyperbolic systems 
of field laws, thus assuming important properties such as finite propagation 
speeds of shock waves and well posedness of the Cauchy problem. The closure 
\ of the system of balance equations is obtained by imposing the entropy princi- 

ple and that of galilean relativity. If we take the components of the mean field 
as independent variables, these two principles are equivalent to some conditions 
on the entropy density and its flux. The method until now used to exploit these 
| conditions, with the macroscopic approach, has not been used up to whatever 

order with respect to thermodynamical equilibrium. This is because it leads to 
several difficulties in calculations. Now these can be overcome by using a new 
method proposed recently by Pennisi and Ruggeri. Here we apply it to the 14 
moments model. We will also show that the 13 moments case can be obtained 



O 



Or 



> 



■ from the present one by using the method of subsystems. 
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1. Introduction 

The 14 moments model was firstly investigated by Kremer [1], up to second 
order with respect to equilibrium; here we want to exploit it up to whatever 
order. The appropriate balance equations for this model reads 

t F + d k F k = 
d t F l + d k F lk = 
d t F lj + d k F^ k = P <ij> 
d t F m + d k F mk = P iU 

where the independent variables are F, F'\ F^ , F M , F vM , which are symmetric 
tensors. See also ref. [2] for further details. The right hand sides of eqs. (JT}i,2 
are zero, such as the trace of that that in eq. dTJ) 3 for the conservation laws 
of mass, momentum and energy. The entropy principle for these equations, by 



using Liu's theorem [3], ensures the existence of parameters called Lagrange 
multipliers, or mean field, such that 

dh = XdF + X,dF l + A, ,'//•"■' + X m dF iU + \ iM dF uu 

d<j, k = XdF k + X,dF kl + A, ,'//'''■' + X m dF kw + X m dF km 

o = \,,l> " + X m P M + X m P m > 0. (2) 

Following the idea exposed in ref. [4], we take the components of the mean 
field as independent variables and define 

h! = XF + X l F i + A, ; /'•'••' + X lU F M + X m F m - h 

<j)' k = XF k + X t F kt + XijF kij + X u F kUl + X lM F ktUl - <\> k . (3) 

By differentiating eqs. ([3]) and using eqs. (H])i,2 we obtain 

dh' = FdX + F l dX,. + / '-'-/A,, + F M dX iU + F iai dX llU 

d<j)' k = F k dX + F kl dX t + F ktj dX l3 + F kai dX m + F kutl dX im . (4) 

In the next section a new methodology recently proposed by Pennisi and Rug- 
geri [5] will be applied (see also [6]) to investigate eqs. ([3]) together to those 
expressing the Galilean Relativity principle, showing that they are equivalent 
to the subsequent conditions (fl"4f . (| 1 5|) and (fT2]) . The last one of these will 
be investigated in section 3 while the other two in section 4. To this end we 
will need the expansion of h! and <f> k up to whatever order with respect to 
equilibrium; it will be introduced also in the next section. 
In section 5 it will be shown how the 13 moments model can be obtained as a 
subsystem of the present one. 

In section 6 we will see that the results of the kinetic approach are a particular 
case of those here found with the macroscopic approach. 
Finally conclusions will be drawn. 

2. The Galilean relativity principle and the entropy principle 

We want now to impose the galilean relativity principle. To this end we recall 
firstly how variables transform with a change of galileanly equivalent frames 
with relative velocity v. For the independent variables, from refs. [1], [2], [7] 
we have 

F = m 
Fi — mi + mVi 
Fij = rriij + 2m^Vj^ + mviVj 

Fui = mm + m n vi + 2m u vi + m^ 2 + 2miviVi + mv 1 v i 
Fan = m uu + Amum + 2m u v 2 + imuVzVi + Arriiviv 2 + mv 4 ; (5) 



here the m... are the tensors corresponding to F... in the second reference frame. 



Moreover we have 

F k = Fv k + m k 

F lk = FiVk + m tk + m k Vi 

Fijk = FijVk + rriijk + 2m k ( i v j ) + m^Vj 

Fuik = F iU v k + m iUk + m ku vi + 2m ku vi + m ki v 2 + 2m k iViVi + m k v 2 Vi 

Fuuk = FmiVk + muiik + 'imkiuVi + 2m ku v 2 + Am kH viVi + Am k iviv 2 + m k v A 



h = h 



hv k + (f> k . (6) 



The first two of these, as the trace of the third and fourth ones, are identi- 
ties, while what remains is the transformation law of the dependent variables. 
Substituting the relations above into eq. ((2])i and defining 

A = A + + X tj ViV.j + \ ipp v 2 v l + Xp Pqq v 4 

Xi — A l + 2X lj v j + 2Xj pp VjVi + XippV 2 + 4:Xppq q v 2 Vi 

Xij = Xij + XhppVhSf + 2XippVj + 2X p pqqV 2 Sf + 4:XppqqViVj 
Xffl Xipp -\- 4:XppqqVi 

Xppqq — Xppqq (7) 



we have 



dh = Xdm + Xidrrii + A,, <////,, + Xuidmm + Xuiidmm ■ (8) 



For eq. ©6 we note that eq. © is the counterpart of eq. in the second 
frame; this allows us to see that eqs. (O are the transformation rules for the 
Lagrange Multipliers. 

Similarly, by substituting eqs. ([6]) in eq. we find 

d(p k = Xdm k + Xidm ki + X i3 dm kij + X M dm kM + X li udm ki iu (9) 

which is the counterpart of eq. ([2])2 in other frame. 
The counterparts of eqs. ([3]) in the second frame are 

h! = toA + rriiXi + ii>, ,\, , + m nU X % u + mmXuu - h 
4> k = 7n k X + m ki Xi + m kij Xij + m ki iiXiu + m k iiuX li u ~ (f) k -, (10) 

differentiating them and using eqs. ©6.7, © and (0 we obtain respectively 

dh' = mdX + rriidXi + mydAy + m in dXui + muudXuu, 
d<f> k = m k dX + m ki dXi + m ki3 dXij + m ki udXiu + muuidXuu. 

Taking their derivatives with respect to the various components of the main 
field we have 

dh' dh! dh' dh! dh' 

m = — — , rrii = m tj = — =— , mm = — » — , mm 



dX dXi dXij dXm dXuu 

d$ k d$ k d$ k d$ k d$ k 

mu = — m ki = — t— , m kij = — r— , m kiU = — — , mum = — ~ — • 
dX dXi dX^ oXui dXuu 

(11) 



Comparing the correspondent terms in the two rows of eq. pip we obtain the 
following compatibility conditions: 

dh' _ d$ k dh' _ d$ k dh' _ d$ k j 

dX k dX ' dX k i dXi ' dX iU dXij 

^=0, -4^ ^=0. (12) 



dXiy dXkm dXm d\u[i 

By substituting h, h, (f> k and <j) k from eqs. Qi, |l0])i. d3j>2 , OH 2 into eqs. 
©6,7, these become 

h' = h', 4 k = 4>' k + h'v k , (13) 
where @i-5 and ([7]) have been used. 

Now, from eqs. (fT3|) we see that h! and cj> k are composite functions of h and 
(f)' k and of eqs. (J7J); but ft,' and depend only on A, A,, Ay , Ay;, A^z and not 
on Vh- In other words, the derivative of h! and k with respect to through 
the above mentioned composite functions, must be zero, i.e. 

dh 1 - 

- — = = mX h + ZrritXih + X ipp (mji^ + 2mjft) + Am hn X ppqq (14) 

= = m k Xh+2m ki X ih +X ipp {m k u8 i + 2m kih ) +4m khU X ppqq + 5%ti (15) 
where (fTT|) and ([7]) have been used. 

The entropy principle and that of material objectivity reduce in imposing eqs. 
(fT4|) . (fl"5]) and (fl"2j) . We want to impose these conditions up to whatever order 
with respect to thermodynamical equilibrium. This is defined, see [8], as the 
state where all the components of the main field, except A and Ajj = j^XuSij, 
amounts to zero. To avoid an excessive quantity of indexes, we will do later the 
expansion with respect to X ppqq . The expansion of the tensor <j> 1 with respect 
to the other variables is 



OO OO OO -. 

\ " \ " \ " J ) iii--ipji---jqkihi—krh r \, \., 



niq'.r 

p=0 q=0 r=o ^ H 



■X jq ii ■ 



X kl h x — -^Xll&kihy j •■' (x kr h r — -^XllSk r h r ) (J ! >! 



With ^■^■■■^ lhl -.. krhr ^ A;;, Xppqq) = 

/ QP+q+rl'i 



dXi ± ■ ■ ■ dXi p dXj 1 u ■ ■ ■ dXj q iidX kl h 1 ■ ■ ■ dX/ lrkr 



(17) 



Now, from the compatibility conditions p2[) 9 . (fl"2")) fi and §F2$ a we see that we 
can exchange the index i respectively with each other index taken from ii ■ ■ ■ i p , 

jli ' ' ' jq and hi ■ ■ • h r or ki • ■ • k r , so (f>p) 1 ,r P ^ 1 J '* 1 1 r r is a symmetric tensor 
with respect to any couple of indexes. Moreover (ftp^q^ 1 3q 1 1 krhr depends 
only on scalars, so that 

f ^ q ;^-^h,- kThT = Q . f p+q+2r+1 is odd 

{^^■■■^^■■■^ = <j,p <q<r (X, Xu,Xp pq q)6 iil ■ ■ ■ S k ^ if p+q+2r+l is even, 

(18) 



so that (j)p^q,r p ^ 1 1 1 r r is known except for a scalar function. 
Similarly, for the tensor h' we have 

oo oo oo 

h '= E E E ^^r^ klhM ^ ■ ■ ■ M*« ■ ■ ■ Ku ■ 

p=0 q=0 r=0 

Afci/u — -^^llSkiln ) ' • ■ ( Xk r h r ~ 7;Xll5k r h r ] (19) 



with h^P 1 "* hlhl " k * h *(\, X u ,X ppqq ) 
( d p+q+r 'h' 



(20) 



V d\i 1 ■ ■ ■ dX ip dXj 1 u ■ ■ ■ dXj udXfahx ■ ■ ■ dXh r k r I 

Taking the derivatives with respect to Xju of the compatibility conditions ([12")^ 
and (|12|> 9 and using (|12|) k we see that we can exchange every index taken from 
JW'jq with each other. Similarly, taking the derivative of eq. (|12[h with 
respect to X rs and using eq. (|12[k we see that we can exchange every index 
taken from ii, - ■ ■ ,i p with each other. Consequently, tfpjq,r P ^ J '« 1 1 fcr ' lr jg a 
symmetric tensor with respect to any couple of indexes; moreover it depends 
only on scalars, so that 

f ^^-i^-krhr = Q j f p+q+2r ig odd 

\^ q p^^--- k ^ = h p . q . r (X, X lh X pPqq )5 1112 ■ ■ ■ 5 k ' K if p+q+2r is even. 

- ...-kH-kh (21) 
In other words, also h p \ q ,r FJl q 1 1 is known except for a scalar function. 

We want to avoid to use eqs. (fT7|) and (|20| in the sequel. To this end we note 

that we can consider 



and 



dXi 1 ■ ■ ■ dXi p dXj 1 u ■ ■ ■ dXj q udXk 1 h 1 ■ • • dXh T k T 

Qp+q+rfik 



dXi 1 ■ ■ ■ dXi p dXj 1 u ■ ■ ■ dXj q iidXk 1 h 1 ■ • • dXh r k r 

depending on X a b as composite functions through A< a b> = (s^Sl — ^S^dabj Xij 

and A;;. With this in mind let us take their derivatives with respect to A a b, 
after that contract them with S a b and calculate the result at equilibrium; we 
find 

7 i-L—i p ji--j q k 1 h 1 —k T h T ab ? _ o_^_uii---i p ji---j q k 1 h 1 ---k r h r / w \ 
n p,q,r+l ° ab — a\ P,q, r v / 

OXu 

and 0*;+^ - j ^ hl - k ^ ab 6 ab = ^^^■^■■■i,^h 1 ...k r h r ^ (23) 

dXu 



An interesting consequence of eq. ([22]) can be observed as follows. 

Let us take the derivative of h! with respect to Ay taking into account that 

Xij = T^XuSij + X < ij > : 



f)%l oo oo oo „, l 

OAij p=0 9=0 r=0 CA ;; 

Afci/i! - -XuSfohi J ■ ■ ■ (^k r h r — g Xu5k r h r ^j % + 

oo oo oc 

EEE -^^-^^-^4, • • • A, A, „ • • • X hU ■ 

p=0 g=0 r=l ^ y ' 

Afci/ii — gAijfJfc^i^ • • ■ yXk r -ih r -i — -^Xu6k r ._ 1 h r - 1 ^J (&hr^kr ~ ~^h r k r ^ 

which, by using eq. (l22|) . becomes 

_p=0 q=0 r=0 



9Ay 



x. x x i i 



\ ^ \ " \ ^ 1 1 1 ii—ipj'i— j q k%hi—k r h r ab x \ \ 



AjiH ' ' ' Aj 9 zz ( A/jj/jj — -XuSkiin ) • ■ • ( Afc r ? lr — -\u5k r h r ) 5 13 + 



3 " ^"W V r T 3 



oo oo oo .. 

r 1 



u ■ 



EEE ^^'^ ' ''' hr h • • • h^u ---K 

p=0 q=Q r=l 

Afei^i — -^Xu5k 1 h 1 ^j • • • ^Afc r _ 1 / lr _ 1 — -\u6k r _ 1 h r _ 1 ^J Sh r k r 8 l] + 

oo oo oo 

+ E E E ^ ft*;^-^-^^ ■ ■ • a,.a ;;; • • ■ hu ■ 

p=0 q=0 r=l 

^A fel h 1 - -XuSk^i^ ■ ■ ■ (^ kr _ lhr _ 1 - -A/ i (5fe r _ 1 / lr _ 1 ^ Sh r §i r ; 

We note that the term in square brackets amounts to zero as can be easily 
proved by substituting r=R+l in the second sum. What remains can be written 

as 

dh' dh' 



dXij 9A < jj> 

where the derivative in the right hand side has been taken without consid- 
ering that the components of A<y> aren't independent because restricted by 
A<ij>^ y = 0. Proceeding similarly with <fi k and using eq. (|2"3")1 we find that 

dcf) k d(j) k 



dXij 9A<y> 

After that, we see that eq. (fT7|) and (|20| become consequences of eqs. (fT6|) and 
(fT9|) so that they can be forgotten. But, instead of them, we have to impose 
eqs. ((22J) and J2SJ) . 

Expliciting eq. (|22|) by means of eq. (|2"Tj) we have 

_ P + g + 2r + 1 dh p ^ r 
W+1 - 3 P + 9 + 2r + 3 e Xu ' (24) 



from which 

hpqr = y p + q + 1 / r V , (25) 

p ' 9 ' r p + q + 2r + l QXlt K ' 

as it can be seen by using the iterative procedure. 
Similarly, expliciting eq. ([23]) by means of eq. (T8)) . we have 

, p + q + 2r + 2d(j> p . q . r 

^ r+1 - 3 p + q + 2r + 4 d X u ' (26) 



from which 

_ or p + q + 2 d r 4>p, qfi 
^ W P + q + 2r + 2 d y a ■ 

that can be proved using the iterative procedure. 
If we introduce the quantities 

{k p ,q = hp^ q fl if p+q is even 
k p ,q = <P P} q,o if P+q is odd, 



(27) 



(28) 



we note that h! and <j> k are known if we know all the terms of the infinity 
matrix k p . ? ; so our aim is to find fc p , g . We have also to impose the compati- 
bility conditions (fT2"|) and the conditions (Lf| and (fTS"]) expressing the Galilean 
relativity principle. Let us begin by investigating the conditions (|12p . 

3. Exploitation of the conditions (1121) 



Now let's impose conditions (fl2|) on our tensors. We notice that equations 
([T2j) 4 k are already satisfied because the tensors 0p,V,'r tp31 3q 1 1 r r are sym- 
metric, so there remains to impose eqs. (fT2"|) i 9 3 5. 



Eq. (p])i. by using (16]), (18]), (20]) and (21]), becomes 
which, for r=0 reads 



hp+l.q.r — (29) 



«p+l,q,0 - -5 (.OUJ 

OA 

and, for the other values of r is consequence of (2"5"]) . (|27p . (3D]). This last one, 
by using (28]), can be written also as 

kp+i^q = — with p+q+1 even. (31) 

In other words, the elements with p+q+1 even of the matrix fc p +i, 9 can be 
expressed in terms of that of the same column but previous row. 

• Let us impose now eq. (|T2]) 9. using eqs. (T6| . (T8]) . (|20| and (21]); we 

obtain 

hp^q^r+l = 4>p+l,q,r (32) 



which, by using eqs. (|25j) and (27| is equivalent to 



and this, by using 



, becomes 
P + <7 



1 dk n 



P 



3 dX u 



with p+q even. 



(34) 



Using (|3ip or (13"4")) we can express all the elements of the matrix k PtQ in terms 
of those in the same column and previous row. Iterating this procedure each 
element can be expressed in terms of the elements in the first row of the matrix. 
In fact joining eqs. (|3"Tj) and (|3"4"|) we obtain 



(k - qg 9+1 

'u.o = 3 2 - 



>,9 



g+2 



3* 



3 £+i^+l 



with p and q even, 
with p and q odd, 
with p even and q odd, 



(35) 



j^rko^q with p odd and q even. 



P+ q+2 OXu^OX^- 



• Finally, let us consider eqs. (TT!?)h 5 . 
they become respectively 



P- 



and 



By using eqs. 
tively into 



and 



p + q 

p + q- 



Using eqs. IP ]) . <[T5 |) . (p0|) and 

(36) 



2r 



-2r 
2r 



-2 
3 



9A fcfc « p + q 

(ITT)) and finally 



2r + 1 



(37) 



>,<?+! 



dkp q 

1 — o — - — 

j + dk 



the above equations transform respec- 
with p+q+1 even (38) 



3 8X 



with p+q+1 odd. 



(39) 



aabb 



each element of the matrix k PtQ can 



In other words with eqs. (f3"5| and 
be written in terms of the element in the same row and previous column. But 
we already know, by eqs. ([35]) . each row of the matrix k P:q in terms of the 
first one; so we have to investigate the compatibility of these two results. By 
substituting eqs. ([35]) into eqs. (138]) and (|39]l we obtain a series of equations 
for the first row of the matrix fc„ „, i.e.. 



ko,q+l — 3-^—k() q 
9+1 d 2 u "_ d h. 
^q+ZdW , q ~ d\ °- q+1 



2+1. 



9+3 d\ aa . 
d d 



q odd, 
q even, 
q even, 
q odd, 



(40) 



and other equations which are consequences of these last ones. Now eqs. (|4"0)) i 
and (|40p ^ give each element fco.g in terms of fco,Oj i- e - 5 



. I- 1 



1 d q 



k , q = 3 2 — — — ; —k 0i0 with q odd. (42) 

q + 2 d\ 2 dX a 2 ahb 

Through these two equations is possible to express a generic element in the 
first row in terms of the first element in the same first row of the matrix. 
Eq. ([40]) 9 a remain to be imposed. The first one of these with q=0 and by use 
of J3TJ reads 

d 2 d d 
9 T^7 fc o,o = — — *o,o, (43) 



dXf, ' dX dX 

aabb 



which is a condition on fc ,o- After that eq. (|4H)> 9 for the other values of q is a 
consequence of eq. (|4"5j) . 

At last, eq. (|4"0)U with use of (|41 [) and (|42p becomes equivalent to its value for 
q=l, i.e., 

Q3 _ d 3 ] 

aabb 



OXldX aabb 

dXdXl 



which is eq. (|43[) differentiated with respect to X aa bi,; so it is sufficient to impose 
eq. 63]). 

We can now substitute eqs. (|4"Tj) and (|42[) into eqs. (|55|) which now become 

fcp.g = i p+g 3 ^" : - fc o.o with p and q even, 

kp. q = 3^~ -r^n . p+q _j d * v +j . jrrfco.o with p and q odd, 

+ _i P q d\ u ^~d\^2-dX aabb — 

kp.q = 3^-5— -T^Tn . p+g-i ^ p" i 5tt^o,o with p even and q odd, 

fcp.a = 3 2 rr~rn ; p+^+i 8 ". p-i i j-fco,o with p odd and q even. 

(44) 

In this way all the elements of the matrix k p>q are determined in terms of fc .o 
which is restricted, until now, only by eq. (|43|) . Another restriction will be 
found in the next section. 

4. Exploitation of the conditions (|14[) and (|15[) 



There remains now to impose eqs. (|14l) and (|T5]) . but we can see that ([14]) is a 
consequence of (1T51) and (fl~2]) . In fact 



• the derivative of fT4")) with respect to Afc is equal to the derivative of (fT"5)) 
with respect to A, thanks to (fTTj) . (H^i, 

• the derivative of (|14p with respect to A^ is exactly the derivative of (jT5J) 
with respect to A&, thanks to pip . (HU)2,i, 



• the derivative of (fTTj) with respect to Afc;; is exactly the derivative of (|T5f 
with respect to A a fc, contracted after derivation by 5 a b, thanks to (|lip . 

• the derivative of (fT4")l with respect to Afc^n is exactly the derivative of 
with respect to A,-«, contracted after derivation by Ski, thanks to (fTTj) . 



Consequently, eq. (fT4"|) needs to be imposed only for A^ = 0, A b = 0, A feH = 
and Xkkll — 0, and in this case it is an identity. So it remains to impose only 



cq. (|15p . To this end it is useful to use the identity 

d r (« chj^\ . d r+1 4>' k 

Ay "777 — "*7 



d\k 1 h 1 ■ ■ ■ dXk r h r \ d\i J d\id\k 1 h 1 • • • dXk r h r 

d r <h' k 



d\f ll d\k 2 h 2 • ' ' d\k r h r ) 



whose proof can be found in the Appendix of ref. [5] and holds also if, in our 
case, (f> fc depends on the further independent variable X aa bb- 
Let us take now the derivative of eq. (|15|) with respect to A^ • • • Aj , X^u ■ ■ • Xj u , 
Xktht ■ ■ ■ Xk r h r - If we calculate it at equilibrium and we use eqs. (JT7J) and ([20)) 
we obtain 

c $ ,i2—ip)kji—j q hiki—h r k r , 2 - ,khii---i p ji---j q hiki---h r k r 

1 r J h 1 k 2 ---h r k r )ki 1 ---i p j 1 ---j q „ ,khi 1 ---i p j 1 ---j q h 1 k 1 ---h r k r 

r jj2—Jo)fc»i— ip/iifci— h T k T ab x . f ,ki 1 ---i p j 1 ---j q hh 1 k 1 ---h r k r 

+ Q°h(j 1 < Pp,q-l, r +l °ab + ^ A aabbWp,q+l , r + 

+ hi)- q ^ h "- j ^ kl --- hrkr S hk = 0. (45) 
To evaluate this condition it will be useful to do the following considerations: 
1) Let V'"' be a symmetric tensor; it is easy to prove that 

fih(ii 1 pi 2 ---i P ji---j q ei---e a k) _ P gh(ii ^12- ■ -i P )ji---j q ei---e s k _|_ 

_fih(ei q pe2"-e s )ii-"i p ji—j<ik _|_ 



p- 


-9 


+ 


s - 


- 1 






'7 






p- 


-'7 


+ 


s H 


- 1 






s 






p- 


-'7 


+ 


s H 


- 1 






1 







p + q + s + 1 



2) Moreover we have 



,j2---j q kif-i p h 1 k 1 ---h r k r ab f - _ , g + P + 2r + 3 x (j 2 ...j a ki 1 ---i p h 1 k 1 ---h r k r ) 

<Pp,q-l.r+l dab ~^ q - 1 ' r+1 q + p + 2r + l 

3) Finally, we can express everything in terms of the scalar h pqr using the 
following relations: 

gj0p-l,g,r = ^p,g,r from eq. (gOJ, 

0p+l,g,j— 1 = hp,q,r, 4>p+l,q,r = /ip.g.r+l from eq. (j32)). 

0p,g-l,r+l = p+^+2r+3 fe P^^ Ir ° m ec l- (03)' 

^P.g+l.r = itlXltloxlj p^r from eq. © . 

All these results allow to rewrite eq. ([45)) as 



= h p . q . r (p + q + 2r+ ij^fe^fa-fcrfcrfc) + |X H< 5^n-^fe^ Pj9ir+1 + 

„ P + g + 2r+ l J . fcMl ... hrfcr ^ * p + q + 2r + 1 dh p ^ r hkil ... hrkr 

+Zq — i To — r~Q n pqr + 4:A aabb — ■ — — — -— t o , 

p + q + 2r + 3 p + q + 2r + 3 dX aabb 



where the notation S eie2 "' e2s — S^ eie2 ■ ■ ■ 5 e2a ~ ie2!,s > has been used; the result is 
equivalent to 

2 - 

= (p + q + 2r + l)h p ^j. + -Xuh p , q , r +i + 
p + q + 2r + l ( - dh p ^ r 

+ i 7T> T^T Zqh vqr + ^Aaabb — 

P + q + 2r + 3 V d\ aabb 

This equation, by using eqs. ([25]) and f28]). becomes 

d r d r+1 d r dk 

0=ip + 3q + 2r + 3)-^ r k Ptq + 2X ll -^ TI k p , q +4X aabb - 



ax r u dx r + l dx r u dx aabb 

with p+q even. We note that if this relation holds until a fixed r taking its 
derivative with respect to Xu we obtain that it holds also with r+1 replacing 
r. Therefore, it suffices to impose this relation for the lower value of r, i.e for 
r=0. In this case it becomes 

d dk 
= {p + 3q + 3)k M + 2X u - r k p , q +4X aabb -^-, (46) 

oXu dXaabb 

with p+q even. 

Let us firstly analyze the case with p and q even. Putting eq. (|35[h into (|46]l 
we have 

= (p + 3q + 3) „ p k . q + 2 An „ - - k 0<g + AX aabb , p , ,q . 

8X 2 dXi dXV dXi OXZdXi dX aabb 



We note that if this relation holds until a fixed p taking its derivative with 
respect to Xu and then with respect to A, we obtain that it holds also with p+2 
replacing p (p must be even). Therefore, it suffices to impose this relation for 
the lower even value of p, i.e for p=0. 
In this case it becomes 

= (3<? + 3)k a . q + 2AH^-fco,, + 4A aa66 ^2-, (47) 

OXu dXaabb 

that is (|4"6"]) calculated in p=0. 

By using eq. (|4"T]) we see that eq. (|47l) becomes 

Qq ^ Qq 8ko 

= (3<7 + 3)— -j — — &o,0 + 2A;; TTT\ — ~ ^0,0 + ^X a abb—^ — — — ~ — — ■ 

d^dX 2 abb dX 2+1 dX 2 abb dX 2 dX 2 abb dX ccgg 

We note that if this relation holds until a fixed q taking its derivative with 
respect to A;; and then with respect to X aabbl we obtain that it holds also with 
q+2 replacing q (q must be even). Therefore, it suffices to impose this relation 
for the lower even order of q, i.e for q=0. In this case it becomes 

= 3k 0fi + 2X u -^k , +4X aabb -J^, (48) 

CM;; dXaabb 



that is (|4*S|) calculated in p=0, q=0. 

There remains the case with p and q odd. We will see that it will give only 



identities. In fact, putting eq. (|35[> 9 into (|4"rj|) , this becomes 



gp 

= (p + 3 9 + 3)— ^ — k , q + 

„« <9 P+1 - 9f <9fc 0g 

d\ H 2 + d\^T- 3\ n 2 dX^ dX aabb 

We note that if this relation holds until a fixed p taking its derivative with 
respect to A;; and then with respect to A, we obtain that it holds also with p+2 
replacing p (p must be odd). Therefore, it suffices to impose this relation for 
the lower odd value of p, i.e p=l. In this case it becomes 

0=(3q + + 2A„- T -^fc , 9 + 4A aa66 °' 9 . (49) 

OX OXoXu OXoXaabb 

This relation, by using eq. (|4"2"j) becomes 

d q d 
= (3g + 4) ^ q+1 - T fc ,o + 

r)X 2 f)\ 2 ^A 

UA U UA aabb 

„t 9 , , ? 99 a 2 /c , 

+ : n _i , : — -m-— F fc 0,0 



9-1,-, .4+1 Q C '"U,u i '"aooo „ 4 _i A g+ i * a- 



bb 11 aabb 



We note that if this relation holds until a fixed q, taking its derivative with 
respect to A;; and then with respect to A aa bf,, we obtain that it holds also with 
q+2 replacing q (q must be odd). Therefore, it suffices to impose this relation 
for the lower odd value of q, i.e q=l. In this case it becomes 

Q2 q3 a <9 3 fcg o 

= 7 ; „ fc ,o + 2A» . „ fc ,o + 4:X ppqq „ ' — , 

dXdXpp qq dXdXudXppqq dXdX ppqq 

which is a consequence of (|4"5)l because it is its second derivative with respect 
to A and X PP qq. In this way, we have seen that the conditions (|14ll and (fT5|) 
give only the restriction (|48| for fco,o and many identities. 
So we have that every element of the matrix k p>q can be expressed 
as function of fco,o and this is restricted only by eqs. (1431) and (|48|) . 
Let us conclude by exploiting these conditions and let us do it by using the 
expansion of fc 0; o around the state with A pp99 = 0, i.e., 

k ,0 = "MX kl)^ppqq- (50) 

Using eq. (|4^|) becomes 

dXl dX 



while eq. (|48|) transforms into 



"~ 1 ^ 1 dk 1 



s=0 s=0 



■KXi 



si - - <—> .s 1 - ^ r..,- m s 



i.e. 




2Aa££n- = fors=0, 
2A„f^+4sfc 5 = fors>l; 



but the relation for s=0 is contained in the other equation, so that they can be 
written as 

dk 

('3 + As)k s + 2X u ^- = Vs>0, 
9Xu 

whose solution is 

k s = A~^%(A). (52) 



This allows to rewrite eq. ([511) as 

^+I=fc 9(3 + 4 S )(5 + 4s). (53) 
dX 4 V 

In this way we have found that fco(A) is an arbitrary single- variable function, 
while the other functions fc s +i(A) are determined by (|53|) . except for a numer- 
able family of constants arising from integration. 

5. The 13 moments model as a subsystem of the 14 moments one 

To verify that the 13 moments case is a subsystem of the 14 moments one 
we will show that the relations obtained in [5] for the scalar functions jo,g are 
satisfied by the value of fco j9 found here but considering X ppqq — 0. Firstly we 
have to rewrite the expressions of ko >q - Substituting eq. (j5Dj) into eq. (j¥T|) we 
have 

_ „| 1 f 1 dik s diX s ppqq _ 

K0,q — 2 — > -7 — — — 

(7 + 1 S! f)\ 2 pl\2 

* s=0 aA U OA ppqq 



1 ^ 1 d2k s q 

7 / ~1 —S(S — 1) • • • IS — - 

^ + 1 ti sl dxl 2 



If we calculate this for X ppqq = 0, only the term for s = | remains, so our 
relations becomes 

g 1 dika 
ko q = 3 2 y~ with q even. 

q+1 

Substituting eq. (|50|) into eq. (J42J) , still making the previous considerations, 
we have 

I d 3 ~2~kg+i 



1 + 2 dX 



ko,q = 3 2 — — q _ 1 2 with q odd. 



7/ 

Now using eqs. (f5"2"|) we obtain 

q 3. - 3 + 3q „. 

35^- (-5) 2 r;(3 + 2<?, 3q + 1)A„ 2 ft a for q even, 

(54) 

3^ ^ (-1)^,7(5 + 2q, 3q + 2)X~^k a±1 for q odd. 



where 77(0, b) = 0(0 - 2)(o - 4) • • ■ (6 + 2)6. 

Comparing this result with the corresponding one for jo,? m [5] (i.e. eqs. (56) 
and (57)), we find that they are the same, except for identifying 

W =(-lY ^T^ 3 + 2 «> 3 1 + l )H W (55) 



and for setting c q = 0. 

It is easy to verify that with I q given by eq. ([55]) , the condition (58) of ref. 
[5] becomes exactly the present eq. (|53|) . except for substituting q=2s+2, and 
viceversa. All the other results of ref. [5], for the 13 moments model, can 
be obtained by substituting X aa bb — in the present ones except for the new 
restriction c q = 0. 

In other words, in ref. [5] the solution was found except for two families of 
constants, one arising from integration of eq. (58) in ref. [5] and another 
constituted by the constants c q appearing in eq. (57). This second family of 
constants doesn't appear if the 13 moments model is obtained as a subsystem 
of the 14 moments one. 

6. The comparison with the kinetic approach 

The solution of our conditions proposed by the kinetic approach, see [2] and 
[9], is 

h! = J F(\ + XiC 1 + Ay-cV + \uic l c 2 + XaabbC A )dcidc2dcz 
'" - F(X + XiC 1 + Xtjc'c' + Xmc'c 2 + X aabb c A )c k dc 1 dc 2 dc 3 , 



(where F is related with the distribution function at equilibrium) , and it is easy 
to see that it satisfies the conditions (|T2|) . (jT4j) , (fTS"]) . We can now see that it is 
a particular case of our general solution. In fact eqs. IT7|) and (|20| now become 



Mx—ifji—jqkihi—krhr 
p.q.r 



c i c ii . . . c *P C h . . . c^c^c^c^ . . . c hr c kr dcidc2dc 3 , 



and it is easy to see that eqs. (|22|) and (|23|) are satisfied. 
Eqs. (HHJ) and flU) hold with 



4tt 



up ' q ' r ~ p + q + 2r + 2 



F^ +r Hx + h u c 2 + A aafcb cV +39+2r+3 ^, 



*w = t ^ ri r F^ +r Hx + \x H c 2 + x^y+^+^dc. 

p + q + 2r + 1 J 3 

The eqs. (|25l) and Q27p are consequences of these. The definitions 
become 



A POO 1 

k p , q = — — / F^ (A + - A H c 2 + X aabb c 4 y+ 3 «+ 2 dci{ p+q is 

p + q+lJo 3 

47T 



p + q + 2j 3 



F(f+«)(A + lx u c 2 + A aafcb c 4 )c p+3<?+3 dcif p+q is odd. 



From these it follows 



f°° 1 
h,o = 4tt / F(X + -X n c 2 + X aabb c 4 )c 2 dc 



and it is not difficult to see that eqs. (|4"3")l and (|¥4")1 are satisfied. 
Proof of eq. 1)48 P needs an integration by parts, as follows 







Fc Z dc+-\ U I F'c^dc + iXaabb 



F'c e dc = 



Fr~dc+ I [ — I r-dr 3 

dc 



Fc 2 dc + Fc 3 



3Fc 2 dc 



which is satisfied because 

lim Fc 3 = 

c — >oc 

We can now see that eq. ([50| holds with 

fc, = Ait 



F^(\+ l -\ n c 2 y s+2 dc, 



of which eq. (|5ip is an easy consequence. 

By using the change of the integration variables c — r)X u 2 , we obtain eq. ([5 
with 



47T 



1 

F W( A +ir, 2 )?7 4s + 2 ^. 



Proof of eq. (|52p needs two integrations by part, as follows 



d 

dX 



k s +i 



4tt 



i 

F( S+2 )(A+ V)r? 4s+6 d?7 
3 



4^+ 1 )(A + i ?7 2 )^ 4 ^ 



6^(4s + 5)F (s+1) (A + -?7 2 )r/ 4s+4 d77 = 



-67r(4s + 5)F (s) (A + - v 2 )-jf s+3 
3 2 



Z" 00 1 
/ -9vr(4s + 5)(4s + 3)F (;5) (A + -r? 2 )?? 



2\ 4s+2 



dr\ 



= ^(4s + 3)(4s + 5)fc s . 
Consequently, the kinetic approach suggest to take 

fc (A) = 47r/ F(X + —r] 2 )r) 2 dr), 



(56) 



which is only a change from our arbitrary function fco(A) to the arbitrary func- 
tion F; moreover it considers only a particular solution of the eqs. (1531) . i.e., 
eq. (|56[) . In this way the numerable family of arbitrary constants arising 
from integration of eq. (|53p doesn't appear in the kinetic approach. Then the 
macroscopic approach here considered is more general than the kinetic one. 



7. Conclusions 



We are very satisfied by the present results because we have found the closure 
of the field equations up to whatever order with respect to equilibrium. This 
was never obtained before in literature. Apparently difficulties in calculation 
becomes very elegant and somehow simpler. Moreover, our model inherits all 
the nice properties of the Extended Thermodynamics, such as to be expressed 
in the form of a symmetric hyperbolic system, to predict finite speeds of propa- 
gation for shock waves, and to guarantee well posedness of the Cauchy problem 
and continuous dependence on the initial data. 
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